Nonequilibrium shot noise spectrum through a quantum dot in the Kondo regime: 
A master equation approach under self-consistent Born approximation 
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We construct a number(7i)-resolved master equation (ME) approach under self-consistent Born 
approximation (SCBA) for noise spectrum calculation. The formulation is essentially non-Markovian 
and incorporates properly the interlay of the multi-tunneling processes and many-body correlations. 
We apply this approach to the challenging nonequilibrium Kondo system and predict a profound 
nonequilibrium Kondo signature in the shot noise spectrum. The proposed n-SCBA-ME scheme goes 
completely beyond the scope of the Born-Markovian master equation approach, in the sense of being 
applicable to the shot noise of transport under small bias voltage, in non-Markovian regime, and 
with strong Coulomb correlations as favorably demonstrated in the nonequilibrium Kondo system. 
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Beyond the average current, shot noise (current fluctu- 
ations) can provide deep insight into the nature of trans- 
port mechanisms [1]. In the past decade, most efforts 
have been devoted to the zero- and low-frequency noise, 
including also the full counting statistics [2]. However, 
even more information is stored in the finite-frequency 
(FF) current noise [3-7]. For instance, the FF noise 
is sensitive to quantum statistics, where a crossover be- 
tween different statistics can be revealed in the frequency 
domain. Also, in the quantum regime, which is defined 
by frequencies higher than the applied voltage or tem- 
perature, the FF noise is a powerful tool to probe the 
characteristic timescales of the system dynamics associ- 
ated with intrinsic excitations and interactions. 

Among the various techniques for shot noise calcula- 
tion (including the counting statistics), the master equa- 
tion approach, particularly its number(n)-resolved ver- 
sion [8-11], might be the most convenient one. However, 
this technique is built largely on the 2nd-order Born- 
Markovian master equation, which limits thus its applica- 
tion only in zero- or low-frequency noise, and under large 
bias voltage. In this work, we will first extend the master 
equation approach beyond these limits, making it highly 
non-Markovian and properly account for the interplay of 
multiple tunneling and many-body correlations. We then 
apply this new approach to the challenging nonequilib- 
rium Kondo system to calculate the FF noise spectrum, 
where a profound Kondo resonance behavior will be re- 
vealed. 

The nonequilibroum Kondo system, with the Anderson 
impurity realized by transport through a small quantum 



dot (QD), has been attracted intensive attention in the 
past two decades [12-22]. Compared to the equilibrium 
Kondo effect, the nonequilibrium is characterized by a fi- 
nite chemical potential difference of the two leads. As 
a result, the peak of the density of states (spectral func- 
tion) splits into two peaks pinned at each chemical poten- 
tial. The two peak structure is difficult to probe directly, 
by the usual dc measurements. Nevertheless, the shot 
noise can be a promising quantity to reveal the nonequi- 
librium Kondo effect, although much less is known about 
it. We notice that results on low- frequency noise mea- 
surements have only appeared very recently [23, 24], 
while so far there are not yet reports on the FF noise 
measurements. A couple of theoretical studies [25-28], 
however, revealed diverse signatures (Kondo anomalies) 
in the FF noise spectra, such as an "upturn" [25] or a 
spectral "dip" [28] appeared at frequencies ±eV/h {V is 
the bias voltage), as well as the Kondo singularity (dis- 
continuous slope) at frequencies ±.2eV/h in Ref. [26], or 
at ±eV/2h in Ref. [28]. Also, it was pointed out in Ref. 
[26] that the minimum (dip) developed at ^eV/h is not 
relevant to the Kondo effect, since in the noninteracting 
case the noise has similar discontinuous slope at -iieV/h 
as well. 

In general we describe a transport setup hy H = 
Hs{aj^,a^) + Hies + H' . Here Hs is the Hamiltonian 
of the central system embedded between two leads, with 
a]Jj (a^) the creation (annihilation) operator of the state 
/i). More specifically, for a small and strongly inter- 
acting quantum dot, with only a single level involved in 
transport to realize an artificial Anderson impurity, we 
have 
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In this model we use /i to label the spin-up ("t") and 
spin-down ("J,") states, and /2 corresponds to the opposite 
spin orientation, is the spin-dependent (single) energy 
level, and U the on-site Coulomb repulsive energy (with 
= ajj^a^ the occupation number operator). The other 
two Hamiltonians, i?res and H', describe the leads and 
their tunnel coupling to the central system. They are 
modeled by, respectively, H^cs = J2a=L.R^k''akbl^k^ak 
and H' = Ea=L,i? E^fc(WA;aji6afe + H.c.) with (bak) 
the creation (annihilation) operator of electron in state 
|fc) of the left (L) and right (R) leads. 

For the study of shot noise, the nonequilibrium Green's 
function based calculation scheme is not efhcient. In 
contrast, an alternative one, say, the particle-number(??)- 
resolved master equation (n-ME) plus the MacDonald's 
formula [29], provides a much more convenient method 
for that purpose. Also, the n-ME is extremely suitable 
for studying the full counting statistics (PCS). To our 
knowledge, the existing n-ME scheme is only precise to 
the Born approximation (BA), i.e., up to the 2nd-order 
expansion of the tunnel Hamiltonian [8-11]. Unfortu- 
nately, however, this type of master equation cannot de- 
scribe the small bias transport, since in this case the mul- 
tiple tunneling process between the system and lead is 
heavily involved. For similar reason, obviously, it can- 
not at all describe the Kondo effect, which is actually a 
consequence of interplay of the multiple tunneling and 
the many-body correlation. Therefore, in order to study 
the shot noise behavior through an interacting QD in 
the Kondo regime, one has to include the effect of higher 
order tunneling process in the master equation. In a re- 
cent work [30] , going beyond the Born approximation, an 
improved scheme under the self-consistent Born approx- 
imation was proposed as follows: 

m = -^Cp{t) - E {[<.4p W] +H-C.}. (2) 

Here we set the Planck constant h = \ and will make fur- 
ther convention in the following for a system of units by 
setting e = fc_B = 1 for the electron charge and the Boltz- 
mann constant. In Eq. (2), we define: ct = -f and — , ct = 
— tr; = aj^, and a~ = a^. Also, the supcroperators in 

Eq. (2) read Cp = [Hs,p\ and 4p (t) = Ea=i.fl -^i'^U^) 

while A:iip{t) = Y.JldTCi"l{t - T){U{t,T)[aZp{T)]}. 

C^Ji, {t — t) is the reservoir correlation function (see Ap- 
pendix A for more details). Very importantly, U{t,T) 
is an effective propagator under the spirit of SCBA, 
which considerably generalizes the iJg-defined free prop- 
agator Q{t,T) = e"*''''-*"'^-' in the 2nd-order Born mas- 
ter equation. The SCBA is implemented by defining 
Pj{t) = U{t, r)[aj^p(r)] (here and in the following we use 
"j" to denote the double indices {v, a) for the sake of 
brevity), and closing Eq. (2) via an equation-of- motion 
(EOM) for this auxiliary object: 

f,,{t) = -iCp,{t) - f dt'4''\t ~ t')pj{t'). (3) 



In this equation the 2nd-order self-energy superoperator, 
Yi^\t — t'), differs from the usual one since it involves 
anticommutators, but not the commutators in the 2nd- 
order master equation (See Appendix A for an explicit 
expression). 

Now we proceed further to construct the particle num- 
ber ("n") resolved SCBA-ME. To be specific, consider 
the reduced system state yO^"^ conditioned on the elec- 
tron number arrived to the right lead, which satisfies 

= -^cp^-^ - { K<pL + -.^'Sr., - 

(4) 

Here A^^^^.^, (0 = /o drCL% (t - r) [pj") (t , r)] , while 
the summation over v makes sense in regard to the 
abbreviation of j = {z^, ct}. In particular, p^^\t,T) 
is the n-dependent version of the quantity pj(t,T) ~ 
h{{t, T)la'Zp{T)], satisfying an EOM according to Eq. (3): 

^ J J J 

(5) 

In this equation we introduced A^'^ (t) = 

E.' /; dt'Clllit - t') [e-^^(^-^')[aZ,pf\t')]]. 

The structure of Eq. (4) follows the same idea of 
constructing the 2nd-order n-resolved master equation 
[8, 11], which is essentially equivalent to the counting- 
field a.pproach [9, 10]. Following [11], we split the Hilbert 
space of the electron reservoirs into a set of subspaces, 
each labeled by n. Then we do the average (trace) over 
each subspace and define the corresponding reduced quan- 
tities as p^^\t) and p'-^\t,T). In Eq. (4), moreover, the 

appearance of p^"^^'' {t, r) is owing to a more tunneling 
event (forward/backword) involved in the process of the 
corresponding terms. These considerations also lead to 
the n-dependence structure of Eq. (5), the EOM of the 
auxiliary quantity pj""* . 

The noise spectrum, Si{uj), is the Fourier trans- 
form of the current correlation function Si{t) = 
{I{t)I{0))ss defined in the steady state. Very con- 
veniently, within the framework of the n-ME, one 
can calculate 5*7(0;) by using the MacDonald's for- 
mula [29]: Si{uj) = 2u}J^dtsm{ujt)f^{n'^{t)), where 

(n^{t)) = = TrE„^^P^"HO: and the n- 

counting starts with the steady state (p). Based on 

Eq. (4), one can express ■^{n^{t)) in terms of -4^j^(i) 
and A^^\At). The former has been introduced in 

RfiN ^ ' 

Eq. (2), needing only to replace p(t) by p. The latter 
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reads A^^^^it) = EJo drC^^U* - r)[N,{t,r)i where 

^jit,T) = Then, the MacDonald's for- 

mula becomes 

Sjic.) = 2c.Im5]Tr{2[^^;)^(^)at - <j^(a.)a,] 
/J 

This result is obtained after Laplace transforming 
^^j^-(i) and A^^l^it). More exphcitly, 

/OO 1 / 

where the Laplace transformation of the steady state 
reads p{uj) = ip/uj, and the propagator lA in frequency 
domain is defined through Eq. (3). On the other hand, 

■^^^\r{^) reads 

In deriving this result, we introduced an additional prop- 
agator through Nj{t,T) = U{t - t)Nj{t), with Njlr) = 
alNir) as the initial condition which is defined by 

^{t) = E„"P^"Ht')- ^(^) and N{uj) can be obtained 
via Laplace transforming the following EOMs. (i) For 
A^(a;), based on the n-SCBA-ME we obtain 

N{t) = --^LN{t) - ^ { [a% A^^^,{t)\ + H.c.} 

/J,fT 

+ E { [^rIp< - ^R^-P^^] + H-C-}- (7) 
('zij For U{lu), from Eq. (5) we have 

/t 

'E{K;kw«M-<k(^K'] (8) 

The self-energy superoperator T,^\t — t') is referred 
to Eq. (A7) in Appendix A for its definition. Similar 

as introduced in Eq. (5), we defined here A^^J-Xt) ~ 

E.' i^dt'd^lit - n {e-^(*-*')K;p,(t')]}. 

For the convenience of application, we would like to 
summarize the solving protocol in a more transparent 
way as follows. First, solve iJ{u!) from Eq. (3) and ob- 
tain p{uj) from Eq. (2); then, extract U{ll!) from Eq. (8) 
and N{io) from Eq. (7). With the help of U{uj), U{uj) 
and N{ijj), one can straightforwardly calculate the noise 
spectrum of Eq. (6). 

Now we return to the Anderson impurity model. Sim- 
ply, there are four states involved in transport: |0), 
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FIG. 1: Shot noise spectrum in the Kondo regime, by vary- 
ing the bias voltage (a), the temperature (b) and the band- 
width {W) of tlie reservoirs (c). We assume = e = fcs = 1 
and use an arbitrary unit of energy in this model simula- 
tion, with parameters as Fl = Ta = F = 0.5, et — ^J. = 
e = —2, and f/ = 6. The bias voltage is defined as usual 
by /i_L = "PR ~ V/2. The Kondo temperature is given by 

= i^/^^^M^^^wp-]' having a value of Tk = 0.144 
for the given parameters. 

t): \ i) and \d), which correspond to the empty, spin- 
up, spin-down and double occupancy states, respectively. 
With respect to these states, the reservoir correlation 
function ci^^ is diagonal, i.e., C^i{t) = 5^uC'^}{t) and 

Fq'^II = fI'^^J^i,. Moreover, using these basis states, we 
can reexpress the electron operator in terms of the projec- 
tion operator form, aj^ = Im)(0| + { — l)^\d){il\, where the 
convention (—1)^ = 1 and (—1)"'' = —1 is assumed. Since 
the shot noise spectrum is defined on the steady-state 
current fluctuations, we need first a solution of the steady 
state {p). In steady state, one can express the key op- 
erator in Eq. (2) as 4^!^ = /_"^ ^ v'^^J {^)U{±^)[ap]. 
Straightforwardly, after some algebra, we obtain [30] 

WH[atp-] = [A+M|m)(0| +<(..)(-ind)(M|] , 
U{"Uj)[a^p] = [\-{ojm{p\ + K-{Lo){-ir\-p){d\] . (9) 

In terms of the matrix elements of p, the specific ex- 
pressions of and are given in Appendix B. 

Substituting A^^p, with the result of Eq. (9), into Eq. (2) 
one can first obtain the steady state. Then, following the 
solving protocol outlined above, the noise spectrum can 
be carried out. 
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In Fig. 1 we display the symmetrized shot noise spec- 
trum in Kendo regime (the numerical results are pre- 
sented with the use of fi. e = fc^ = 1). First of all, we 
notice a remarkable dip behavior (Kondo signature) in 
the noise spectrum at the frequencies w = i.V/2, as par- 
ticularly demonstrated in Fig. 1(a) by altering the volt- 
ages. We attribute this behavior to the emergence of the 
Kondo resonance levels (KRLs) induced at the Fermi sur- 
faces, i.e., at = V/2 and = —V/2. In steady state 
transport, it is well known that the KRLs are clearly re- 
flected in the spectral function, i.e., the effective density 
of states (DOS) of the Anderson impurity. In terms of 
the master equation (sec Appendix B), the KRLs struc- 
ture is hidden in the self-energy terms, which characterize 
the tunneling process and define the transport current. 
Similarly, the noise spectrum is essentially affected, par- 
ticularly in the Kondo regime, by the self-energy process 
in frequency domain based on the same master equation. 
This explains the emergence of the spectral dip appearing 
at the same KRLs (i.e., at w = ±.V/2). 

However, we would like to remark that the dip behav- 
ior is also a consequence of highly non-Markovian treat- 
ment of the current correlations. We have checked that, 
using the quantum-jump technique [32] or the quantum 
regression theorem [33], this behavior cannot be recov- 
ered, even the evolution during (0,t) is treated as non- 
Markovian based on Eq. (2). The point is that the defini- 
tion of the current in the correlation function (/(t)/(0)), 
in the non-Markovian case, cannot be independent of the 
propagation during the time interval {0,t), because of the 
non-Markovian memory effect. In contrast, based on the 
n-SCBA-ME, the MacDonald formula correctly accounts 
for the correlation between the current and the memory 
effect during (0, t), by employing the number(n)-counting 
technique. 

Alternatively, as a heuristic picture, one may imagine 
to include the KRLs as basis states in propagating p{t), 
which is implied in the current correlation function. In 
usual case, when the level spacing is larger than its broad- 
ening, the diagonal elements of the density matrix decou- 
ple to the evolution of the off-diagonal elements. How- 
ever, in the Kondo system, the diagonal and off-diagonal 
elements arc coupled to each other, through the com- 
plicated self-energy processes. This feature would bring 
the coherence evolution described by the off-diagonal ele- 
ments, with characteristic energies of the KRLs and their 
difference, into the diagonal elements which contribute 
directly to the the second current measurement in the 
correlation function (/(t)/(0)). Then, one may expect 
three coherence energies, ±V/2 and V, to participate in 
the noise spectrum. Indeed, the dip emerged in Fig. 1 
reveals the co/ierence-induccd oscillation at the frequen- 
cies ±y/2, while the other one at the higher frequency V 
(observed in Ref. [28] in the case of infinite U) is smeared 
in our finite U system by the rising noise with frequency. 

Physically, the current fluctuation spectrum corre- 
sponds to electron transfer between the dot and leads, 
accompanied by the energy (w) absorption/emission of 



detection. Therefore, as the frequency (w) matches the 
energy difference between the dot level and the Fermi 
surface of the lead, certain "singularity" associated with 
the Fermi function at the Fermi surface is expected to 
emerge in the spectrum. This is reflected in Fig. 1(a) 
by the staircase behavior. This "singularity", however, 
has been smoothed by the finite temperature effect (see 
Fig. 1(b) for further illustration). In Fig. 1(c) we dis- 
play the bandwidth effect. For finite (narrow) band- 
width, the spectrum would diminish at high frequencies 
(when much higher than the bandwidth), since in this 
case the electron transfer channel associated with the w- 
emission/absorption is switched off. In the low frequency 
regime, on the other hand, we find that the narrowing 
bandwidth would shift the Kondo dip to lower frequency. 
This feature indicates that the Kondo peak pinned at 
the chemical potential is only a result in the wide band 
limit. For finite (especially narrow) bandwidths, it may 
need further work to determine the location of the Kondo 
peaks. 

To summarize, we have applied a new shot noise 
scheme to the nonequilibrium Kondo system, for finite 
U and arbitrary bandwidths. The scheme is based on 
a generalized number (n)-resolved master equation under 
self-consistent Born approximation, which considerably 
goes beyond the scope of the usual 2nd-order Born mas- 
ter equation. This treatment allows us to predict a pro- 
found nonequilibrium Kondo signature in shot noise at 
frequencies associated with the chemical potentials. We 
anticipate a wide range of applications of the proposed 
approach to shot noise studies, as well as future work to 
clarify the diverse Kondo signatures in noise spectrum 
[25-28]. 
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Appendix A: Some Particulars in the SCBA-ME 
Approach 

1. Reservoir Spectral Density Function 



The imaginary part, through the dispersive relation, is 
associated with the real one as 

A(±) (uj) =V r — ^ r(^) (w) 



in Eq. (2) read 



The key operators 
E«=L,H'4i';U0, and A^^l{t) = EJodrC'aMt 
t) {U{t, T)[a" p{t)]} . cj^J^{t — T) are the correlation func- 
tions of the reservoir electrons (in local equilibrium), be- 
ing defined as 



Re 



^ f /3(c.-,0 
2 27r 

UJ - fla\ 



2 27r y ^ Wa 



(A6) 



C^il(i-r) = (/(^)W/f;(r)) 



B- 



(Al) 



Here, 



.Mvv - and f^'^>{t)^ ^ U,it), 

resulting from rewriting the tunneling Hamilto- 
nian H' = Ec.=L,i? EA'fclWfc^^^A'fe + H.c.) = 
Ea=L,RE,, (aj^fa^ + H.C.), by introducing = 
Efc^apfc^apfe- The time dependence of the operators in 

ci^l/ — r) originates from using the interaction picture 
with respect to the reservoir Hamiltonian, while the 
average {■ ■ ■)b is over the reservoir states. Moreover, we 
introduce the Fourier transform of Cc^Ju (t — r) through 



dio 
2^' 



(A2) 



where V stands for the principle value and is the 

digamma function. 



2. Anomalous Self-Energy Superoperator 

The central idea of the SCBA-ME scheme is replacing 
the free propagator in the 2nd-order master equation, 
Q{t,T) — e-»^(*— by an effective one, U{t,T) under 
the SCBA spirit. By introducing Pj{t) ~ U{t,T)[a1 p{t)], 
we obtain Eq. (3), the EOM of this auxiliary object. In 

Eq. (3), the 2nd-order self-energy superoperator, Y}2^\t — 
t'), is worth receiving some special attention. As labeled 
by the superscript "(A)", an anticommutator, instead of 
the usual commutator, is involved there. That is, the 
self-energy superoperator has the following form: 



Accordingly, we have ri|it(a;) 
and T^aJuioJ 

27r tc^fet* — efc) is the spectral density function 

of the reservoir (a), ri'a\^) denotes the Fermi func- 
tion na{uj), and ni \uj) = 1 — naiuj) is introduced for 
brevity. Alternatively, we may introduce as well the 
Laplace transform of C^Ju{t — t), denoting by C^Ji,{uj), 

which is related with ri^ly(i^) through the well known 
dispersive relation: 



(w) - T^,,{co)ni-\oj), where r„,,,H = / di'S^^^^lt - = E } + {4' <^ } 



where A^^^ is defined as 



e-*^(*-*'n<:ft(t')]}- 



2tt cj ± cj' 



iO+ 



Tt^iico'). (A3) 



In this work, for the reservoir spectral density function, 
we assume a Lorentzian form as 



(A4) 



In some sense, this assumption corresponds to a half- 
occupied band for each lead, which peaks the Lorentzian 
center at the chemical potential pa ■ Wa characterizes the 
bandwidth of the ath lead. Obviously, the usual constant 
spectral density function is recovered from Eq. (A4) in the 
limit Wa — > oo, yielding Taiiu{i^) = ^a^iv Correspond- 
ing to the above Lorentzian spectral density function, 
straightforwardly, we obtain 



We remark that the anticommutative brackets appeared 
in Eq. ( A7) indicate that the propagation of pj {t) 
does not satisfy the usual 2nd-order master equation. 
This actually violates the so-called quantum regression 
theorem. 



3. Steady-State Current 

Similar to the usual 2nd-order master equation ap- 
proach, the current through the ath lead reads 



Icit) 



2e 



E Re {Tr[^W (t)a^ - ^^(t)at ] } 



(A8) 



citl(^) = i [rw (T^) + zAW (T^) 



(A5) 



Moreover, the steady state together with its associated 
current can be obtained easily as follows. Consider 
the integral Jq dT[- ■ ■]p{t) in A'^^pit). Since physically, 

the correlation function C^u{t — r) in the integrand is 
nonzero only on finite timescale, we can replace p(r) in 
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the integrand by the steady state p, in the long time limit 
(t oo). After this replacement, we obtain 



A. 



(±) 

a^^p 



E 



dw 
2^ 



rW(c.)i^(±^)[a±p]. (A9) 



Then, substituting this result into Eq. (2), we can 
straightforwardly solve for p and calculate the steady 
state current. 

We would like to mention that, remarkably, for nonin- 
teracting system, the steady state current given by this 
SCBA-ME scheme coincides precisely with the nonequi- 
librium Green's function approach, both giving the exact 
result [30]. Notice also that, by contrast, the Born mas- 
ter equation is applicable only to sequential tunneling 
transport, being valid only in large bias limit. 



Appendix B: Steady State Solution of the Anderson 
Impurity Model 

In Eq. (9), associated with the steady state solution of 
the Anderson impurity model, we have 



(Bi) 



and r^(^). Now, if is proportional to by a 
constant, the steady state current can be recast to 
the Landauer-Biittikcr type, in terms of an integra- 
tion of tunneling coefficient over the incident energies, 
I = T^^!-o. t {^l{oj) - T'i^)- The tunnel- 

ing coefficient, very compactly, is given by 7"(a;) = 
Tr{rLrR(rL + rfl)-iRe[<^(w)]}, where <p(w) = + 
ip2 (w) . For the Anderson impurity system in nonequilib- 
rium, we find 



+ 



^[^^,^H-4+)M](l-n^) 
»(i - 



(B4) 

where = + pdd, and 1 - = p^^ + poo- This 
result, precisely, coincides with that given by the EOM 
technique of the nonequilibrium Green's function [31]. 
One can check that, as discussed in detail in Ref. [31], 
this solution contains the nonequilibrium Kondo effect. 



Here we introduced n^^(a;) 
and n5-J(w) 



UJ 



energy SomI"^) is given by 

dio' 



Som(^) Theself- 



2tt oj -oj' + iQ+ ' 



(B2) 



while I]J(cli) by 



duj' 



27r UJ 
duj' 



(^') 



iO+ 



271 UJ ^ Ed 



^0^ 



(B3) 



With the above results, as outlined after Eq. (9), one 
is able to carry out the steady state solution p. 
Based on it, to obtain further the current, we first 
introduce (y9i^^(w) = Tr [a^pii,(a;)] and ip2fj.ui^^) = 
Tr[a^iP2iy(w)] , where pi,j{uj) and p2u{^) are calculated 
using Eq. (3), with an initial condition of piv{^) = pa], 
and p2i/(0) ~ o\,P- To simplify notations, we denote 
the various matrices in boldface form: ipi{uj), ip2{i^) 
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